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Abstract 

Data collected around =91 GeV by the OPAL experiment at the LEP e + e~ collider 
are used to study the mechanism of baryon formation. As the signature, the fraction of 
X~ hyperons whose baryon number is compensated by the production of a £~,A or S - 
antihyperon is determined. The method relies entirely on quantum number correlations of 
the baryons, and not rapidity correlations, making it more model independent than previous 
studies. Within the context of the JETSET implementation of the string hadronization 
model, the diquark baryon production model without the popcorn mechanism is strongly 
disfavored with a significance of 3.8 standard deviations including systematic uncertainties. 
It is shown that previous studies of the popcorn mechanism with AA and p7rp correlations 
are not conclusive, if parameter uncertainties are considered. 
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1 Introduction 



The formation of baryons within a jet of hadrons has proved difficult to model and is still not 
well understood. While the shape of momentum spectra can be derived from QCD with the 
modified leading logarithmic approximation together with parton-hadron duality [1], more 
complex observables like correlations have not been derived from first principles. 

Several physical models like the thermodynamic [2], cluster fragmentation [3] or string 
fragmentation [4, 5] models, have been developed to describe baryon production. Of these, 
the most successful is string fragmentation, based on the creation of diquark-antidiquark pairs 
from the vacuum as illustrated in figure 1. In a production chain of particles along the string, 
a baryon and an antibaryon can be produced in immediate succession (figs. la,b) or else one 
or more intermediate mesons can be produced between them (figs. lc-e). This production of 
intermediate mesons, referred to as the popcorn effect [6, 7], is included as an option in the 
Monte Carlo event generators JETSET and PYTHIA [8, 9] and can be steered with a free 
parameter. 

Past experimental investigations of the popcorn effect made use of rapidity ordering of 
hadrons in the fragmentation chain. Intermediate mesons modify the rapidity difference be- 
tween associated baryons and antibaryons. Rapidity correlations between A A pairs produced 
in Z° decays have been studied by several LEP collaborations and the conclusion was that best 
agreement between the experiments and the JETSET Monte Carlo model [8] was obtained 
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Figure 1: Baryon production in the diquark model without (a,b) and with (c-e) the popcorn 
effect. 

with the popcorn effect included [10]- [12]. 

A contradictory result was reported by the DELPHI collaboration [13]. Their measure- 
ment is based on triple correlations between a proton, an antiproton and a charged pion close 
in rapidity. Because the popcorn effect enhances the pion density in the rapidity interval 
between the proton and antiproton, the minimum rapidity difference between a pion and a 
proton was compared for the particle orderings p7rp and 7rpp. The measurement indicated 
that the rapidity rank correlations could be reproduced without the popcorn effect. The con- 
tribution of events with popcorn produced mesons was reported to be less than 15% at 90% 
confidence level. Insufficient modeling of the fragmentation dynamics could not, however, be 
excluded [13]. 

In this analysis, the popcorn mechanism is investigated in a different way, by tagging 
rare baryons and measuring the quantum numbers of correlated antibaryons. By not relying 
on rapidity differences, we obtain results that are more model-independent than previous 
studies. The data were collected with the OPAL experiment at the LEP e + e~ collider at 
CERN. Especially suitable for our purposes is the S~ hyperon. If a E~ is produced, its 
baryon number and strangeness are compensated either by an antinucleon and a kaon, as 
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illustrated in figure la, or by an associated antihyperon. The case of A or E° production 
without the popcorn effect is shown in figure lb. The same graph, with the bottom-most 
u-quark replaced by a d- or s-quark, describes associated S~ and H~ production. Associated 
production of a £+ and H° antihyperon with a S~ can only occur through the popcorn 
mechanism (figs. ld,e). This makes and correlations ideal tools to study the 

popcorn effect. Unfortunately, the production rate of S~ hyperons and the probabilities for 
simultaneous reconstruction of £~ and S+ or E~ and H° particles are too small to make 
such an analysis feasible, given the available data statistics. The equivalent analysis cannot 
be performed with tagged A hyperons because in this case S+ and H° antihyperons can be 
produced without the popcorn effect. Alternatively one can measure the fraction 

F n = F v-jF + F v-Z + F v-?F (1) 

of £~ hyperons accompanied by a S~ (F s - xp)' a ^ ( F n- a) or a "~ ( F z~ These 
correlations can occur in the popcorn model but are more likely in the diquark model (fig. lb) 
and their rate is thus a sensitive measure of the baryon production mechanism. 

The exact definition of the three correlations F s _ ^ in the sum (1) needs to account for the 
possibility that an event may contain more than one £~ hyperon, other additional hyperons, 
or more than one antihyperon. More generally, one can consider an arbitrary particle k and 
its antiparticle k. If the number of S~ hyperons in an event is larger than 1, all combinations 
of E~ k and X~ k pairs are counted. Denoting the total rates of £~-antiparticle and S~- 
particle pairs by i? s _ ^ and -Rs-,fc> respectively, and the total S~ rate by the fraction 

F s _ ^ can be written as 

FE __ = ^i_^ (2) 

This definition implies that the antiparticles k in an event enter n^n s - times. The E~ £~ 
pairs are counted n^- (n s - — 1) times. In the data analysis, this multiple counting is not 
an issue because the number of reconstructed S~ hyperons per event is very small. Baryon 
number conservation ensures that Sfe-^E-fc = 1' ^ the sum extends over all antibaryons, 
including the antihyperons and antinucleons. It is the understanding throughout this paper 
that the charge conjugated channels are included. Thus, the rates R^- ^ contain all like-sign 
pairs £~ k and £~ k, the rates ^ contain all unlike-sign pairs E~ k and S~ k, and £~ 
antihyperons are included in . 

It is evident from figure lc that kaons created through the popcorn effect reduce Fjj-. 
The PYTHIA Monte Carlo program without the popcorn effect, tuned to reproduce the 
observed baryon rates and momentum spectra, predicts F-jj « 0.9, as will be shown later. 
The other extreme is a model in which baryon number and strangeness are compensated 
statistically, i.e. from conservation laws alone. The ratio of weakly decaying hyperon to total 
baryon production in Z° decays is approximately 0.22 [14]. Because there must be at least one 
antibaryon in the rest of the event, if the E~ is detected, and the tagging biases the number of 
strange valence antiquarks, the lower bound will be a bit larger: F-jj > 0.22. A more rigorous 
calculation can be performed with the thermodynamic model for particle production. The 
advanced version of the model, based on the microcanonical ensemble, yields F-jj = 0.23 for 
an initial system without strangeness [15], close to the simple estimate of the lower bound. 
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2 Event selection 



2.1 Event topologies 

In any experiment with a sufficiently large tracking device S~ hyperons can be identified 
by track kinks from X~ — > n-zr" decays. At LEP energies, the efficiency is much less than 
100%, because the decay vertex lies often outside the fiducial volume for its reconstruction. 
In principle, correlated decays S~ — > n7r + and E~ — ► A7r + could be reconstructed using track 
kinks, too, but the efficiency for the exclusive reconstruction of £~ S _ or S~ H~ pairs is 
small. 

In this work, only the £~ hyperons were reconstructed exclusively. Two signatures for 
correlated antihyperons were used: 

1. A hyperons were reconstructed by analyzing their so-called V° topology from the decay 
A — > p7r + . The decay vertex and the flight direction of the A candidates allow the 
A impact parameters do with respect to the beam line to be computed. Direct A 
production, including the contribution of decays from S° hyperons, is characterized by 
low impact parameters do, while large impact parameters indicate a preceding weak 
decay and are a signature for S decays. 

2. Charged pions with significant impact parameters are a signature for weak decays of 
arbitrary antihyperons. An inclusive sample of tracks with large do values, consistent 
with a pion interpretation, was selected. Throughout this paper, this data set is referred 
to as the sample of displaced tracks. This sample has a large background, and different 
antihyperon species contribute with different weights, because the number of decay 
pions per antihyperon is 1 for £ - , 0.64 for A and 1.64 for chain decays of S _ . 

The correlated S~A candidate sample gives the numbers of true £~ A and £~ H~ pairs. 
The fraction F^p can be extracted from the displaced track sample by a weighted subtraction 
of the £~ A and S~ H~ contributions. 

2.2 Experiment and data sets 

All data taken by the OPAL experiment [16] in the Z° energy region during the years 1991 to 
2000 were analyzed to measure the correlation. The OPAL experiment had nearly complete 
solid angle coverage and excellent hermeticity. The innermost part of the central tracking 
detector was a high-resolution silicon microvertex detector, which immediately surrounded 
the beam-pipe [17]. It was followed by a high-precision vertex drift chamber, a large- volume 
jet chamber [18], and z-chambers, all in a uniform 0.435 T axial magnetic field. In this work, 
the outer detector parts as well as the forward detector system were needed for triggering 
and identification of multihadronic events only. The criteria for multihadronic event selection 
have been described elsewhere [19]. 
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The present analysis is entirely based on the central tracking system [16, 18]. For candi- 
dates to be accepted, all central wire chambers and the microvertex detector were required 
to be fully operational. The data sample for this analysis consists of 3.97 million events. 

The identification of strange particles is based on earlier work [20, 21]. The innermost 
sense wires of the jet chamber had a distance of 25.5 cm from the beam spot, the wire-to- wire 
distance was 1 cm and there were 159 sensitive radial layers. The requirement of the pattern 
recognition program was the existence of at least 12 hits in the jet chamber. This makes it 
possible to identify S~ hyperons with decay lengths larger than 36 cm. 

The z coordinate, along an axis parallel to the electron beam, was measured with a 
precision of 700 fim with the stereo wires of the vertex chamber and 100 to 350 fim with the 
z chambers. For tracks leaving the drift chamber at the side cones, the z coordinates of the 
exit points can be computed from the radius of the last wire with a hit. At the beam spot 
a constraint can be set using the bunch length of the beam. Inside the jet chamber, the z 
coordinates were measured with the charge division method with a resolution of 6 cm. This 
is one limiting factor for the kinematical reconstruction of S~ hyperons to be discussed later. 

The quality of the impact parameter measurement is directly connected to the detector 
resolutions in the (r, ip) plane perpendicular to the beam axis. These resolutions were 5 fim 
to 10 fim for the microvertex detector [17], 55 /jm for the vertex chamber and, on average, 
135 fim for the jet chamber [16]. 

To study the detector response, the selection was applied also to Monte Carlo samples 
used before at OPAL. These were generated with the JETSET7.3 and JETSET 7.4 programs, 
followed by a full detector simulation [22]. The steering parameters for the generator are 
given in [23, 24]. The subset of parameters relevant for this analysis is discussed in section 4.2 
and the numerical values are given in table 8 (appendix). In total, the Monte Carlo samples 
consisted of 4.65 million multihadronic Z° decays. 

2.3 £ Selection 

The production cross section of S~ hyperons has previously been measured by OPAL [20]. 
Here the selection criteria to find track pairs forming a decay vertex were slightly changed 
relative to [20] to improve the sensitivity of the analysis to correlations. The applied cuts 
are summarized in table 1. They accommodate pattern recognition tolerances and define a 
fiducial jet chamber volume to guarantee minimal hit numbers for both tracks and to remove 
background from the end plates. The ratio h/pi in the last line of table 1 is proportional to 
the decay time in the rest frame of the decaying particle. High values are rejected to reduce 
the substantial background from kaon decays. 
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variable 


condition 


transverse momenta of the decaying 
and the daughter particle 


PTi > 0.15 GeV/c, pT2 > 0.10 GeV/c 


total momentum of the decaying particle 


Pi > 2 GeV/c 


impact parameter of the primary track 


do < 5 cm 


endpoint radius of the primary track 


te\ < 170 cm 


maximum gap between the primary and 
the secondary track in the (r, if) plane 


d\2 < 10 cm 


number of secondary tracks 


1 


particle charges 


equal 


arc distance between the intersection point 
and the end points of tracks 1,2 in (r, if) 


Hi > — 7 cm, fi2 < +7 cm 


radius of the intersection point 


r V tx > 35 cm 


distance of the intersection point 
from the end plate in z direction 


Azytx > 40 cm 


X 2 for agreement of the two tracks in the 
(r, f) plane, using (3 x 3) error matrix 


X 2 > 2000 


ratio of the track length to the 
momentum of the decaying particle 


ii/pi < 15 cm/(GeV/c) 



Table 1: Selection of track pairs forming a decay vertex (for a more detailed description of 
the variables, see [20]). 



Monte Carlo studies show that the following processes have to be considered as sources 
for S~ candidates, others being negligible [20]: 

1. E~ — > nir~; 

2. S+^nvr-; 

3. S - -► Att~; 

4. — ► 7r~+neutrals; K~ — > /i"+neutrals or K — > e _ +neutrals; 

5. secondary reactions in the detector material and fake background. The last sample 
consists mainly of scattered particles. 

To determine the relative contributions of these processes to the observed data sample, an 
unfolding procedure very similar to that used in [20] was used. The invariant mass mj- and 
the decay angle 6* of the pion candidate in the rest frame of the hypothetical £~ particle were 
computed, assuming that the unseen neutral particle is a neutron. The angle 0* is defined 
with respect to the flight direction of the £~ candidate at the decay vertex. 

Figure 2 shows scatter plots of these variables for the data and three Monte Carlo sources 
of events with kinked tracks, namely events with H~ or K in the final state. The 
background 5, not shown in figure 2, is a smooth function of both observables and is largest 
in the backward direction cos 9* = — 1. 
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Five two-dimensional regions, denoted by a, ...e, were introduced in the (m s - , cos 9*) plane 
to enrich dedicated kink sources (see figure 2). The bin number of a X - candidate was used 
as an observable in the unfolding procedure to be described in section 3. 

In total, 16790 candidates were found in the two-dimensional plane in the mass range 
from threshold to 1.5 GeV. In the Monte Carlo sample 18754 kinks were identified. 




Figure 2: Definition of two-dimensional bins to disentangle X - hyperons from background. 
The reconstructed mass is plotted versus the pion emission angle in the hypothetical center- 
of-mass system. Top left: data. Other plots: Monte Carlo simulations for three particle 
classes as indicated. The labeling of the two-dimensional bins is shown in the bottom left 
plot. The preferential bins for the kink sources are: a+d for S~ and S+, b for b,c and d 
for K~ and e for background. 
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2.4 Selection of pions from weak decays 



Two pre-cuts were applied to select tracks with a good reconstruction quality: the transverse 
momentum with respect to the beam had to be larger than 0.15 GeV/c and the track angle 
at the beam spot relative to the beam direction was restricted to the region | cos6*| < 0.80. 

To remove charged particles from charm or bottom decays, the impact parameter d§^ 
with respect to the primary vertex was required to be larger than 0.2 cm. This cut is the 
essential condition to define the sample of displaced tracks. A lower value would be sufficient 
but does not improve the accuracy of the correlation analysis. 

Pions were enriched using the specific energy loss measurement of the central drift cham- 
ber [26]. A weight w^E/dxi^) was defined as the probability that the energy loss dE/dx of a 
pion deviates from the median value < dE/dx > (ir) by more than the measured difference 
from the median value. The applied condition was w^E/dxi^) > 0.02 and the number of hits 
contributing to this measurement had to be at least 20. 

Finally, the angle between the momenta of the pion candidate and the £~ candidate at 
the primary vertex was required to be less than 90 degrees. This hemisphere cut is motivated 
by the fact that it rarely happens that baryon number and strangeness are compensated by an 
antihyperon in the opposite event hemisphere. The cut reduces the combinatorial background 
by a factor 2. 

In total, 9965 correlated like-sign £~7r~ and 11951 unlike-sign £~ ir + pair candidates 
were selected with these cuts. The corresponding results for the Monte Carlo sample are 
10769 and 13818, respectively. If the total number of Monte Carlo track kinks is scaled to 
the observation, the number of like-sign pairs in the data is well reproduced by the Monte 
Carlo, the difference being (+3. 4 ± 1.4)%. The observed unlike-sign minus like-sign difference, 
however, deviates from the prediction by (—28 ±7)%. This deficit, already visible at raw data 
level, indicates that the Monte Carlo sample contains too many correlated antihyperons and 
is the basis for the final result of this paper. 

2.5 A Selection 

The selection cuts to find A decays in the central drift chamber have been described in [21]. 
The preselection cuts were relaxed from those of [21]. In this analysis all candidates with 
reconstructed masses between the p7r threshold and 1.20 GeV/c 2 were accepted. This larger 
mass window was needed to study the non-A background. 

Two criteria were added. If the reconstructed A flight path points back to the (r, ip) 
position of the hypothetical £ decay kink within 2 degrees, it was assumed that the kink 
originated from a H~ — ► Air~ decay and the A candidate was dropped. The cut removed 
approximately 3/4 of the like sign H~ A (S~ A)-pairs and reduced the self-correlation of H~ 
particles with their own decay A's accordingly. 

Secondly, the hemisphere cut applied to the decay pions was also applied to the A's. The 
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angle between the flight directions of S and A at the primary vertex was required to be less 
than 90 degrees. 

For the correlation analysis, the mass window was reduced to a ±10 MeV wide interval 
around the true A mass. Totals of 276 S~A (£ _ A) and 604 S~ A (£ _ A) pair candidates 
passed all selection cuts. In the Monte Carlo sample 284 like-sign and 721 unlike-sign pairs 
were found. The observed unlike-sign minus like-sign difference is smaller than the Monte 
Carlo prediction by (16 ± 11)%, if normalized to the number of observed kinks. 



3 Correlation analysis 
3.1 Unfolding of kink sources 

From the Monte Carlo sample one gets, for every kink source i and every two-dimensional 
(m s - , cos #*) bin j, the number of accepted events K^ MC \j). The populations of the bins 
j =a,...e in figure 2 are sensitive to the invariant mass and 9* resolutions, which in turn depend 
on the z resolution of the central drift chamber. The z coordinate, however, is not well modeled 
in the Monte Carlo simulation. Because the track end points in the jet chamber are well known 
by other measurements, this mismodeling at the decay vertex gives the dominant contribution 
to the systematic error of the E~ rate [20]. To correct the Monte Carlo program for it, the z 
components of all S momenta were modified according to p z , n ew = Pz,rec + c- (p z ,rec — Pz,true), 
where p z ,true is the true momentum from the MC generator and p z ,rec the reconstructed 
momentum. The constant c is one common factor to be determined in the analysis. 

The total number of kinks K(j) expected in bin j is given by 



where A^ ata and Nmc are the total number of multihadronic data and Monte Carlo events, 
respectively. Incorrect Monte Carlo rates are corrected for by the five scaling factors £j. If 
they are known, the true production rates per multihadronic event, Ri, can be computed for 
all sources i, for example 



Contrary to [20] the factors £j were treated as momentum independent, because the analy- 
sis of [20] had shown that the modeling of the spectral shape was satisfactory. The contri- 
bution from £+ hyperons to K(j) is less than 20% of that of the E~ particles [20]. The 
ratio of the genuine production rates, approximately 1 due to isospin symmetry, was fixed 
to the Monte Carlo prediction, so that £|pr = £s- . Four scaling factors £j were thus left for 
adjustment. They were computed with a x 2 -^ of the measured kink rates to equation (3) 
for a given value of c. 




(3) 



*e- = &--4 mc) with 4^ = 



(4) 



j 
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It was then checked whether the Monte Carlo simulation reproduces the reconstructed 
mass and cos 9* distributions. Differing from our previous analysis [20], the proportionality 
factor c was chosen to get the lowest x 2 -sum for both distributions. The best overall agreement 
was found for c = 1.35 with a one-sigma interval ranging from 1.32 to 1.45. The final fit result 
is shown in figs. 3 and 4. The agreement is excellent. The error of c is not included in the 
statistical errors of the £j. It is treated separately as a contribution to the systematic error. 




1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 



Reconstructed E" mass (GeV/c ) 



Figure 3: Measured £~ mass spectrum. Points: data. Histograms: results of the fit, ordered 
according to the sources for track kinks. The plotted errors are the statistical errors of the 
data. 
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The X - + X~ production rate was found to be R^- = 0.073 ± 0.004. According to [20] a 
systematic error of ±0.009 has to be added. Within the total error, this result is consistent 
with our published value R%- = 0.083 ± 0.011 and also with the world average R%- = 
0.082 ±0.007 [14]. The difference is mainly due to the modified treatment of the z-resolution. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos(0 ) 



Figure 4: Cosine of the center-of-mass pion emission angle. Points: data. Histograms: results 
of the fit, ordered according to the sources for track kinks. The plotted errors are the statistical 
errors of the data. 
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3.2 Correlated particle sources 

The correlated rates of kink-V° or kink-displaced-track pairs were measured as a function 
of the bin number j and the impact parameter do of the correlated particle relative to the 
beam line. The variable do contains information on the lifetime of the parent particle of the 
correlated particle in case of a preceding weak decay, and is needed to disentangle correlated 
A's from S's decaying into A's. 

In the Monte Carlo simulation, nine sources have to be considered for the displaced tracks 
and correlated V°'s: 

1. TF->n7r+; 

2. £+ — > n-7r~ or p7r°; 

3. S~ — > Att + — ► p7r + -7r + ; 

4. H° -> Avr° -> pvr+vr ; 

5. A — > p7r + , including A antihyperons from S° decays, but without the contributions from 
H decays; 

6. charged particles from K + decays; 

7. charged particles or V° configurations from K° decays; 

8. non-K° background of the V° topology; 

9. displaced tracks from secondary interactions, mainly scattering. 

All sources, except for the eighth one, contribute to the sample of displaced tracks. For V°-like 
events only the sources 3,4,5,7 and 8 are relevant. The sample of displaced tracks is enriched 
in pions but contains also a small fraction of protons and leptons. These contaminations are 
included in the Monte Carlo rates and are classified according to the above scheme. 

The combination of five kink sources with nine correlated particle sources leads to a total 
number of 45 classes of particle pairs. In addition, one has to distinguish between like-sign 
and unlike-sign pairs. Equation (3) can be generalized to the like-sign pair rate D like (j,d ) 
and the difference of the unlike-sign pair rate D unhke (j , do) and D hke (j,do), 

nliket- j \ _ ST-J 1 ) N data n {MC,like) , . , x (t .s 

D U, «o) - 2^ %k ■ ' D i,k U>«oJ , (5J 

i,k 

D unlike {jj do) _ D like {jj ^ = £ ^«) . ^ata . ^MC, unlike) ^ ^ _ fl (JC W ^ ^ (g) 

i,k 

The indices at the Monte Carlo rates f)^ C ' bke ) anc [ jj{mc, unlike) g p ec jfy ^ e ^j n ^ source % 
and the correlated particle source k, respectively. The parameters rjf^ k and rjf^ are 90 scaling 
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kink source 


correlated particle 


E- 


S+ 




E° 


A 


K+ 


£- 


0.33 


0.014 


0.15 


0.01 


0.28 


0.04 


S+ 


0.01 


0.32 


0.01 


0.14 


0.28 


0.024 




0.34 


0.026 


0.15 


0.029 


0.36 


0.45 


K 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


0.03 


0.50 



Table 2: Fractions F^ k as predicted by the PYTHIA 6.1 generator with the parame- 
ters from ref [24]. 



factors. In principle, they have to be extracted with a combined fit of the j and do dependent 
distributions of kink-track pairs and kink-V° pairs to equations (5) and (6). 

The Monte Carlo generator predicts the generic correlations F:~ , defined with equa- 
tion (2), generalized to arbitrary sources i and k. Table 2 gives the results for hyperon- 
antiparticle correlations at generator level for the most recent parameter set used by the 
OPAL experiment [24]. 

If the scaling factors rlf^ are known, the experimental results for the correlations can be 
computed 

Here, the scaling factors £j from the kink fit (section 3.1) enter. 

The factors rjfj, are only needed to parameterize the statistical errors for the fit using 
equation (6). With very few exceptions, they are close to one and no physical result is 
extracted from them. 



3.3 Evaluation of the correlation matrix 

In spite of the large number of scaling factors rjf^ , a reliable fit of the data can be obtained. 

This can be seen from the integrals of the differences (D^ C ' unhke \j, do) — D^ C ' hke \j, do)) 
over the variables m£-,cos#* and do, listed in tables 3 and 4. The entries in the tables 
are normalized to the total differences, using rffl = 1, and thus show the relative importance 
of the terms. The nine most significant correlations are given in the tables, the individual 
contributions of all other sources being less than 2%. In total, the entries in table 3 for 
displaced tracks account for 96% of the total rate difference. The contributions to the V° in 
table 4 add up to 104%, the excess being compensated by a small amount of E + A correlations 
with the opposite sign. 

The scaling factors for the four largest correlations in tables 3 and 4, r/°2 T , r/°2 = — 

and rj^l -~ are determined with the fit. Three of them are needed to compute Fjj. 
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kink source 


displaced track source 


E- 




A 


K+ 


background 


£- 


(27 ±1)% 


(16 ±1)% 


(10 ± 1)% 




(17 ±4)% 




(12 ±1)% 


(6 ± 1)% 


(5 ± 1)% 




K 








(3±1)% 


kink background 











Table 3: Largest contributions to the correlations between kinks and displaced tracks. 



kink source 


V u source 




A 


V° background 


E- 


(20 ± 2)% 


(43 ± 3)% 


(14 ±3.5)% 




(6 ±1.5)% 


(21 ± 2)% 



Table 4: Largest contributions to the correlations between kinks and V° candidates. 



The data statistics do not allow to fit more than four parameters. The remaining 77's were 
thus fixed by symmetry considerations, isospin invariance or Monte Carlo studies. Systematic 
errors were assigned to them, if necessary. In the following, a few examples are described, 
preferentially the correlations in tables 3,4. 

For symmetry reasons, one has 

( a ) ( a ) to\ 

= ^- s^ • ( 8 ) 

There is no model independent prediction for F~_ — . This parameter was set to the original 
Monte Carlo result. It follows then from equation (7) that 

*4-,iF = k-- ( 9 ) 

It has been checked for various acceptable Monte Carlo generator tunings, described in the 
next section, that this procedure is valid within 20%. 

Non-negligible parts of the correlation are introduced by the correlated particle back- 
grounds 8 and 9; they are given as sums over the kink sources in tables 3 and 4. The origin 
of this effect is charge conservation in the events. Since the displaced-track-background con- 
sists mainly of scattered particles, it reflects the original particle charges. The kink selection 
introduces a charge bias for the rest of the event, which is visible in the ensemble of remain- 
ing charged particles on a statistical basis. This is true not only for scattered particles, but 
also for asymmetric fake V° candidates. Kinks and correlated particles are assigned to each 
other by chance. The condition T^ rack = i^fyo = was introduced, assuming a correct 
modeling of the background sources 8 and 9. This is justified because the equivalent relation 
for the like-sign background rates was confirmed with a fit to the data with the parame- 

terizations rjfjrack = ^iCtrack and rjf^ = £,iCv° > the results being (track = 1-01 ± 0.04 and 
Cyo = 1.03 ± 0.07, respectively. 

The remaining entry in table 3 is the charged kaon-kaon correlation. Its contribution is 
small because kaon decay inside the jet chamber is unlikely due to the long kaon life time. Since 
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the charged and neutral kaon production rates are almost equal and the hyperon rates are 
much smaller than the kaon rates, the fraction F K - K + is close to 0.5 so that rj^l_ K+ = Ck-- 

In total, the contribution of all other sources to the observed correlation is smaller than 
the statistical error of the final result. Nevertheless, all sources were investigated in detail to 
minimize the systematic error. 

A correlation potentially dangerous for the fit is the S~K + correlation. The corresponding 
asymmetry in table 2 is small, but actually it is the difference of much larger components. 
Strangeness conservation requires 

F E",S^ + F E",E+ + 2 • F E-,H- + 2 • F E-,^ + F E-,A + 2 • ^-, K+ ~ 1 • (10) 

The parameters F^_ — and ^ appear with the weight 2 because the H-particles carry 

two units of strangeness. The contribution from hyperons is negligible. The factor 2 in 
front of F^- K + accounts for the K° contribution, which cannot be measured. Monte Carlo 
simulations with different model parameters showed that relation (10) is fulfilled within 3%. 
Together with the condition for baryon number conservation, 

^E-.E^ + F E-,E+ + F S-^ + F E-,H0 + F E-,A + F T.~ ,p + F Y.- ,n ~ 1, i 11 ) 

equation (10) gives the relation 

F Z-,K+ ~ 2 " ( F V-,p + F T,~,n ~ F ^-^nF ~ F T,-,SP) ' ( 12 ) 

which shows the presence of large compensating terms. In the analysis, F%- k+ was therefore 
computed with the sum rule (10), taking the small correlations F^_ ^qr and F^_ ^ from the 
Monte Carlo generator. Similarly, F s - K + can be constrained by the equivalent equation for 
the E~, 

F ~-,X= + F H-,S+ + 2 • + 2 • F ~-0* + F H-,A + 2 ■ F Z~,K+ ~ 2 . (13) 

The correlations of the S + hyperons are related to those of the T,~ hyperon by the isospin 
symmetry, for instance F E+ ^ = F E _ No model independent predictions exist for F K _ 
and the popcorn specific correlations F^_ ^ and F^_ which were determined with the 
Monte Carlo event sample. 

The above relations allow either to replace rj^ by the four factors to be fitted, or to fix it 
and its contribution to the right hand side of equation (6). A simultaneous binned x 2 fit, using 
equation (6), was performed for the do distributions of the correlated-track and A candidates 
in the five (m^- , cos #*)-regions. In the A case, the reconstructed mass was restricted to the 
narrow interval given in section 2.5. In parallel, a fit of the like sign pairs to equation (5) was 
performed to determine Ctrack = vf track/ & an ^ three normalization factors for the hyperon-A, 
K~-hyperon and (E~ or E~) — (S~ or £ _ ) correlations. The statistical errors of the pair 
rates, needed to compute x 2 > were computed with the equations (5) and (6); they depend 
on the result of the fit. The fit was therefore done iteratively, setting the unknown r/- factors 
to unity for the first iteration. The statistical errors had do dependent fluctuations due to 
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the limited Monte Carlo statistics. These were reduced with a one-dimensional smoothing 
algorithm described in [25]. The asymmetry fit resulted in a x 2 value of 144 for 146 degrees 
of freedom. The coefficients obtained are correlated, the largest correlation coefficient being 
«-0.7 between — and r/ a ^ — . 

'E-E- 'E-E- 



3.4 Experimental results and systematic errors 

The adjusted do distributions are shown in figs. 5 and 6. The histograms give the contributions 
of the correlated particle sources; all kink sources and two-dimensional bins j =a,...e are 
combined. The errors of the data points are statistical and the corresponding errors of the 
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Figure 5: Impact parameters of the displaced tracks, correlated to S~ candidates. The 
difference between unlike-sign and like-sign combinations is shown. Points: data. Histograms: 
results of the fit, ordered according to the sources for displaced tracks. The plotted errors are 
the statistical errors of the data. 
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Monte Carlo histograms are not shown. Fig. 7 gives the rate differences between the unlike- 
sign and like-sign kink-V° pairs as a function of the reconstructed A mass. Both distributions 
are very well described. 



g 100 
3 



80 



60 - 



40 



20 







-20 



OPAL data 



H _|_ M 



0.001 0.01 



D V° background 




d of correlated V (cm) 



Figure 6: Impact parameters of the A candidates, correlated to S~ candidates. The difference 
between unlike-sign and like-sign combinations is shown. Points: data. Histograms: results 
of the ht, ordered according to the V° sources. The plotted errors are the statistical errors of 
the data. 
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Figure 7: Reconstructed masses of the A candidates, correlated to £~ candidates. The 
difference between unlike-sign and like-sign combinations is shown. Points: data. Histograms: 
results of the fit, ordered according to the V° sources. The plotted errors are the statistical 
errors of the data. 
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correlation 




F — 

E-,E- 








fit result 


0.162 ±0.100 


0.054 ± 0.047 


0.234 ±0.057 


0.449 ± 0.091 


systematic error 




mismodeling of correlated 


+0.009 


±0.003 


±0.014 


±0.026 


particle momenta 


±0.005 


±0.002 


±0.007 


±0.013 


z-resolution at 
kink vertex 


_i_n nQ7 


_i_n noQ 


_i_n ri/i k 


_i_n no/i 

±U.UZ4 


dE/dx calibration 


±0.010 


±0.005 


±0.004 


±0.008 


hemisphere cut 


±0.005 


±0.002 


±0.007 


±0.013 


number of hits per track 


±0.007 


±0.004 


±0.002 


±0.006 


cos 9 distribution of tracks 


±0.012 


±0.001 


±0.001 


±0.012 


A detection efficiency 


±0.006 


±0.004 


±0.004 


±0.002 


charge asymmetry 
of detection efficiencies 


±0.002 


±0.001 


±0.001 


±0.002 


uncertainty of 
popcorn channels 


±0.008 


±0.005 


±0.004 


±0.007 


uncertainty of 
correlation 


±0.003 


±0.001 


±0.001 


±0.003 


uncertainty of 
K~A correlation 


±0.001 


±0.001 


±0.003 


±0.002 


final result 


0.17 ±0.11 


0.057 ±0.056 


0.25 ±0.08 


0.48 ±0.10 



Table 5: Experimental results and systematic errors. 



The fit results are listed in the first row of table 5. The errors are statistical and include 
the data and Monte Carlo contributions. Internally, the fit gets the A part of the correlation 
essentially from the correlated A sample, which sets also bounds on F^_ — . The 
correlation is computed from the displaced-track sample as a difference. Neither F s _ |p nor 
— differ from zero in a statistical way, in contrast to the overall sum Fjj. The smaller 
error of the sum is due to the strong anti-correlation between F^_^= and The final 

results, including the systematic errors, are given in the last row of table 5. 

The systematic errors, listed in table 5, will be discussed in the following. 

Mismodeling of correlated particle momenta. The momentum spectra of the correlated 
pion and A candidates are not well reproduced by the Monte Carlo program. The simu- 
lated momenta have to be scaled downwards by as much as 20%. The effect exists both 
for the like sign and the unlike-sign track and V° candidates. As a consequence, the 
measured fractions F { ^ are systematically too small. To find a correction, the detection 
efficiencies for the correlated particles were extracted from the Monte Carlo sample and 
the shifted momentum spectrum was folded with the efficiency function. This leads to 
an upwards correction of the measured correlations. A global correction of (6 ±3)% was 
applied. This correction, together with its error, is much larger than the uncertainties 
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due to miscalibrations of the particle momenta or the impact parameters. Therefore, 
no additional errors were assigned to the momentum and do selection cuts. 

z-resolution at kink vertex. In the ratio of a correlated rate to the single S~ rate the 
overall detection efficiency for £~ hyperons cancels. However, an uncertainty of the 
number of observed S~ hyperons arises from the unfolding procedure. The dominating 
error source is the mismodeling of the z-resolution. The correction factor c, introduced 
above, was conservatively varied between 1.25 and 1.5. For the individual correlations 
.F s _ £ the obtained shifts were found to be non-parabolic functions of c. The maximal 
shifts of the results F s _ -r and Fjj were taken as uncertainties. 

dE/dx calibration. Any antiproton impurity in the displaced track sample reduces the 
asymmetry from direct or indirect A decays. For a pure pion sample and a perfect 
calibration, the frequency distribution of the weight WdE/dxi^)-, defined in section 2.4, 
should not depend on w^E/dxi 71 )- A superimposed peak at w^E/dx^) = is expected 
due to non-pions. The shape of this peak does not perfectly agree with the Monte Carlo 
prediction. A systematic error was assigned to the corresponding mismodeling of the 
antiproton rejection efficiency. 

Hemisphere cut. In principle, the results are corrected automatically for the hemisphere 
cut. No significant correlations were observed in the dropped hemisphere. A systematic 
error would appear, if the fragmentation model were incorrect. A 3% error was assigned 
to all correlations, based on the error of the number of tracks in the omitted hemisphere. 

Other systematic errors The number of hits per track is well modeled and introduces an 
error of 1% for the number of displaced tracks. Also the acceptance cut for the 9 angle 
plays a minor role only. From the difference between the angular distributions of data 
and Monte Carlo events an efficiency error of 3% was estimated. Systematic errors due 
to the cuts for A selection were already discussed in [21]. An uncertainty of ±3.3% 
was taken from that paper as fully correlated error for the A event sample. Differences 
Aej between the detection efficiencies for particles and antiparticles would lead to a 
spurious asymmetry, if both the kinks i and the correlated particles k are affected. 
The effect is proportional to Ae^- • Ae^. The efficiency differences were extracted 
from the observed particle and antiparticle rates and the upper limits for the spurious 
asymmetries in table 5 were obtained. Furthermore, there are uncertainties due to the 
model dependence of the correlations which had to be subtracted. One half of the Monte 
Carlo baryon antibaryon pairs are accompanied by a popcorn meson. The r]^ -factors 
for the popcorn specific correlations and were varied between and 

2 and the fit repeated. The K~A correlation was varied by 50% and an uncertainty of 
20% was assigned to the correlation, fixing the small residual contribution from 

H~ self correlations. As already mentioned, the definition of Fjj introduces a small 
amount of double counting, both in the data and in the Monte Carlo simulation. This 
problem was studied with a toy Monte Carlo program and the only effects found were 
negligible corrections to the statistical errors of the fit. 

The correlated kink-track and kink-V° pairs should be concentrated in the S~ enriched 
bins a and d in the (m s - , cos 6*) plane shown in figure 2. As a cross check, a simpler analysis 
was performed, where the analysis was restricted to region a. All background asymmetries 
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and also the H~A correlation were subtracted as predicted by the Monte Carlo program. The 
final result including reevaluated systematic errors, Fjj = 0.472 ± 0.155, is fully consistent 
with the main result and demonstrates the absence of anomalies in the (m s -,cos6>*) plane. 

4 Comparison with fragmentation models 

Before comparing the result for Fjj to the predictions of models, we adjusted the models, 
incorporated into the PYTHIA Monte Carlo event generator, to describe a set of observables in 
Z° decays. The optimization began with the PYTHIA steering parameters given in [24]. These 
were slightly modified to reproduce the newest experimental information on the baryonic 
sector, including data on AA correlations. The tuned Monte Carlo was then used to predict 
F-jj. We also used the models to study p-7rp correlations more extensively than previously, as 
described below. 

4.1 Observables 

The input observables for the tuning fall into four categories: 

1. Eight baryon multiplicities in multihadronic Z° decay: p, A ++ , A, £+ + £ _ , £°, E~, 
S* + + and The production rates were taken from the compilation of the 
particle data group [14] . 

2. Proton and A momentum spectra. Baryon spectra have approximately Gaussian shapes 
if parameterized in terms of the variable £ = ln(l/x p ) [1]. We use the mean values and 
variances of the proton and A spectra as the observables. Deviations from the normal 
distribution are known, the true maxima being somewhat higher than the result of the 
fit [27, 31]. However, the Gaussian fit is an easy way to compare different experiments. 

The tabulated data [27]-[31] were fitted to Gaussian functions in the interval 1.2 < 
£ < 4.2. The situation for the protons is not satisfactory, the fitted maxima of the £ 
distributions varying from 2.79 [28] to 3.08 [31]. The values i peak = 2.80 ± 0.07 and 
a = 1.11 ± 0.06 were obtained with a combined fit using all LEP and SLD data. The 
errors were not taken from the fit but conservatively estimated from the systematic 
differences between the £ spectra of the experiments. It should be noted that both Cpeak 
and a depend on the fit range due to deviations of the ^-distribution from the Gaussian 
shape. The agreement between experiments is better for A production, leading to values 
e, peak = 2.62 ± 0.04 and a = 1.21 ± 0.04. 

3. AA correlations. The correlation was parameterized by two observables, the first one 
being the rate excess of A A pairs over A A + A A pairs per event, Y^T r '. The other 
observable is the mean rapidity difference Ay A ^- between the A and A after subtraction 
of like-sign pairs. The number of correlated A A pairs per event was taken from [12]. The 
mean AA rapidity difference Ay was computed from the data of [12], as the truncated 
mean for Ay < 3.0. This cut was introduced to suppress the contribution of A's from 
the opposite event hemisphere. 
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4. Nine meson production rates per multihadronic Z° decay: 7r + + tt , it , K + + K , Kg, 
p + + p~, p°, if, K*+(892) +K*-(892) and K*°(892). These rates, taken from [14], were 
included in the tuning to protect the meson generation against parameter modifications 
steering the baryonic sector. 

4.2 Model Parameters 

The popcorn mechanism is incorporated in PYTHIA in two ways [9]. The first simple ver- 
sion was originally introduced in JETSET. Baryon production is controlled by the following 
parameters: 

1. the suppression of diquark-antidiquark production relative to quark- antiquark produc- 
tion, P(qq)/P(q) = PARJ(l); 

2. the suppression of ss production relative to uu production, P(s)/P(u) = PARJ(2); a 
tuning of this parameter was necessary, because the strange meson rates had to be 
readjusted; 

3. a double ratio involving diquarks containing s quarks, (P(«s)/P(ud))/(P(s)/P(u)) = 
PARJ(3); 

4. the suppression factor for spin 1 diquarks, (l/3)P(udi)/P(udo) = PARJ(4); 

5. the popcorn parameter, which determines the relative occurrences of the baryon-meson- 
antibaryon and baryon-antibaryon configurations, PAR J (5); 

6. an extra suppression for having an s s pair in a baryon-meson-antibaryon configuration, 
PARJ(6); 

7. an extra suppression for having a strange meson in a baryon-meson-antibaryon config- 
uration, PARJ(7); 

8. a parameter which enters the exponent of the Lund symmetric fragmentation function 
for diquarks, PARJ(45); this parameter has an impact on the rapidity difference in 
baryon-antibaryon correlations. 

In the advanced popcorn scheme, a universal equation for tunneling from the vac- 
uum is applied to the generation of new partons and an arbitrary number of mesons can be 
created between a baryon and an antibaryon. The tunneling probability is proportional to 
exp(— [3 q • Mj_), where j3 q is a flavor dependent model parameter and M± the transverse mass 
of the created object. Only the first two and the last parameter of the above list are used in 
this scheme. There are three new parameters, two of them related to the tunneling formula: 

9. the tunneling coefficient for u-quarks, (3 U = PARJ(8); 

10. 5f3 = f5 s -f3 u = PARJ(9) 

11. an extra suppression factor for spin 3/2 baryons = PARJ(18). 
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popcorn parameter 
PARJ(5) 




0. 


0.5 


1. 


5. 




observable 


data 


MCI 


MC2 


MC3 


MC4 


MC5 
advanced 
popcorn 


x 2 




56 


59 


59 


63 


99 


j\jcorr 
AA 

^AA 


0.0612 ± 0.0034 
0.71 ±0.04 


0.066 
0.66 


0.060 
0.69 


0.066 
0.67 


0.058 
0.75 


0.081 
0.57 


(this work) 
— (this work) 
F s _ (this work) 


0.17 ±0.11 
0.057 ±0.056 
0.25 ±0.08 


0.39 
0.18 
0.30 


0.34 
0.16 
0.28 


0.30 
0.15 
0.28 


0.19 
0.09 
0.27 


0.20 
0.08 
0.28 


F-jj (this work) 


0.48 ±0.10 


0.87 


0.79 


0.73 


0.55 


0.56 



Table 6: Comparison with fragmentation models. Fit quality x 2 , AA correlations and the X~ 
antihyperon correlations as defined in the text. The errors of the simulation are smaller than 
the last digit shown. 

4.3 Simulation Results 

For a given PYTHIA parameter set, the measured observables from section 4.1 were compared 
with the simulation results and a \ 2 was computed. It cannot be expected from a fragmen- 
tation model that all its predictions are correct to better than a few percent. To reduce the 
contributions of very accurately measured observables, the errors to compute \ 2 were thus 
taken to be at least 2.5%, which represents the characteristic level of agreement between the 
data and MC. Low \ 2 values were searched for with the method described in the appendix, 
for fixed values of the popcorn parameter. Many tunes of the generator have almost the same 
quality. Some PYTHIA parameter sets and the predicted baryon production rates are given 
in the appendix. Table 6 shows the x 2 's, the AA correlations and the S~ -antihyperon corre- 
lations, as a function of the popcorn parameter, where the simple popcorn model is denoted 
MCI to MC4. 

The variable \ 2 is an indicator for the quality of the baryon modeling. For the simple 
popcorn model, the lowest \ 2 value found was 56 for the 23 observables. This means that 
the Monte Carlo generator describes, on average, the observables roughly within 2 times the 
experimental errors or 5%, whichever is larger. 

The most important result of the simulation is that the overall quality of the description 
of the observables does not depend strongly on the popcorn parameter in the simple popcorn 
model. The AA correlation parameters are always reproduced within two standard deviations, 
whether the popcorn effect is switched on or off. The Monte Carlo parameter space examined 
here is larger than that in earlier studies. In view of the overall uncertainty it is not possible 
to reach a definite conclusion about the the popcorn effect by using AA correlations. 

As shown in the bottom line of table 6, the measured fraction Fjj is consistent both with 
the predictions of the original popcorn model with a large popcorn parameter (MC4) and the 
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optimized advanced popcorn model (MC5), but is smaller than the Monte Carlo prediction 
for zero popcorn effect. It was investigated whether models without the popcorn mechanism 
could be found that reproduce the observables of section 4.1 and give Fjj <C 0.9. The relevant 
parameters of the PYTHIA generator from section 4.2, including the diquark fragmentation 
function, were randomly varied as described in the appendix. In these studies, F-jj values less 
than 0.86 were not obtained. The experimental result deviates from this lower limit by 3.8 
standard deviations. 

The advanced popcorn model has fewer parameters available for tuning and provides a 
significantly worse description of data as seen from the larger x 2 value. The larger x 2 f° r this 
model arises to a large extent from two well known facts. The absolute number of A A pairs 
is too large and the distribution of the A A rapidity differences is too narrow in comparison 
to the observation [12]. On the other hand, the average AA rapidity difference in the simple 
popcorn model with a very large popcorn parameter is too large, so that a combination of 
the two models might possibly describe the AA correlation well. 

To complete the comparison of correlations with results from Monte Carlo generators, the 
DELPHI p-7rp [13] correlation was also investigated. This was done for rapidity ordered p7rp 
and 7rpp or pjwr particle configurations inside event hemispheres. The selection cuts and the 
definition of the minimum rapidity gap Ay m i n between a selected pion and the next proton 
were taken from [13]. The discriminating variable of [13] is the ratio of intensities 

R{Aymm) = ^(pvrpj + ^vrpp + ppvr) (U) 

at the rapidity difference Ay m i n . A strong dependence on the popcorn effect had been seen 
by [13] at large values of Ay m i n . For three bins in the range 0.625 < Ay m i n < 1, the observed 
distribution agreed with a subsample of Monte Carlo events without the popcorn mechanism 
and disagreed with a disjunct subsample, containing popcorn mesons, by more than five 
standard deviations, averaged over the three bins. 

These results could be reproduced with the simulations described here: The R distribu- 
tions obtained with the parameter set MCI, but without detector corrections, agreed with 
DELPHI'S observation within 1.9 standard deviations, while there was disagreement between 
the data and Monte Carlo study MC2 by 5.2 standard deviations, averaged over the same 
three Ay m i n bins. Variations of the fragmentation model had not been studied in [13]. An 
increase of the fragmentation parameter PARJ(45) to unity, fixing the other parameters of 
the simulation MC2, reduces the difference between the data and the model prediction to 2.4 
standard deviations. Furthermore, it was found that the advanced popcorn model with the 
parameter set MC5 results in almost the same function R(Ay m i n ) as the model MCI without 
the popcorn effect. 

The modification of the parameter PAR(45) in the simulation MC2 increased the x 2 value 
in table 6 from 59 to 70 due to a shift of Ay^ to a value below the observation, without 
degrading the description of the other observables of section 4.1. In summary, these results 
indicate a high sensitivity of the rapidity correlations to the fragmentation dynamics. 
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popcorn parameter 
PARJ(5) 




0. 


0.5 


1. 


5. 
















MC5 


observable 


data 


MCI 


MC2 


MC3 


MC4 


advanced 
popcorn 


protons 


1.046 ± 0.026 


1.03 


1.03 


1.02 


1.09 


1.01 


A++ 


0.087 ± 0.033 


0.098 


0.118 


0.122 


0.122 


0.127 


A 


0.388 ± 0.009 


0.369 


0.362 


0.366 


0.388 


0.388 


£+ + £- 


0.181 ±0.018 


0.129 


0.131 


0.129 


0.133 


0.141 


S° 


0.076 ± 0.010 


0.069 


0.070 


0.069 


0.070 


0.073 




0.0258 ± 0.0009 


0.029 


0.029 


0.028 


0.029 


0.024 


E*(1385)+- 


0.046 ± 0.004 


0.043 


0.048 


0.052 


0.051 


0.054 




0.0016 ± 0.0003 


0.0005 


0.0004 


0.0004 


0.0004 


0.0004 



Table 7: Baryon rates per multihadronic event. The errors of the simulation are smaller than 
the last digit shown. 



5 Conclusions 

Our investigations indicate that the fragmentation models have not yet reached a state where 
they can quantitatively describe the AA correlations, the pp-zr correlations and the S~ antihy- 
peron correlation Fjj simultaneously. Neither the AA nor the ppzr correlations can provide a 
clear conclusion about the popcorn effect because they can both be described to an acceptable 
level by models either with and without the popcorn mechanism. The PYTHIA generator 
without the popcorn effect can reproduce both rapidity correlations simultaneously within 
two standard deviations. Both correlations are also in acceptable agreement with model 
predictions including the popcorn effect. However, we were unable to find a variant of the 
popcorn model that could simultaneously describe both types of correlations. 

In this work, the mechanism of baryon formation was studied by counting S _ hyperons 
and correlated antihyperons from hadronic Z°decays. The result Fjj = 0.48 ± 0.10, based on 
data taken by the OPAL experiment at LEP, favors a large popcorn parameter in the simple 
popcorn model and is also consistent with the advanced popcorn model. The fragmentation 
dynamics play no role here because the rapidity is not used in the analysis. Correlated particle 
momenta play an indirect role only, because they influence the detection efficiencies, but the 
final result contains a correction and a systematic error for mismodeling. 

Trivial correlations between hyperons and antihyperons based on baryon number and 
strangeness conservation, as predicted by the thermodynamic model, always exist. Due to 
the limited data statistics and the insensitivity to the dynamics, our result gives only a weak 
indication for non-trivial quark correlations between baryon-antibaryon pairs. Statistical 
models, only constrained by conservation laws, differ at most by 2.6 standard deviations. 

The result for Fjj deviates from the lower limit of simulations without the popcorn effect 
by 3.8 standard deviations and thus demonstrates the need for the the popcorn effect in order 
to reproduce baryon correlations within the diquark fragmentation model, where a baryon 
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MP1 


MP 9 


MP3 




advanced 


parameter 














tune 










popcorn 


PARJ(5) 


0.5 


0. 


0.5 


1. 


5. 


- 


PARJ(l) 


0.085 


0.092 


0.100 


0.101 


0.117 


0.226 


PARJ(2) 


0.31 


0.33 


0.33 


0.34 


0.34 


0.31 


PARJ(3) 


0.45 


0.32 


0.41 


0.40 


0.48 




PAR J (4) 


0.025 


0.055 


0.029 


0.021 


0.007 




PARJ(6) 


0.5 




0.33 


0.56 


0.35 




PARJ(7) 


0.5 




0.21 


0.20 


0.28 




PARJ(8) 












1.00 


PARJ(9) 












2.07 


PARJ(18) 












0.18 


PARJ(45) 


0.5 


0.16 


0.22 


0.52 


0.33 


0.27 



Table 8: Parameter sets used for the simulations in tables 6 and 7. 



MC parameter 


PARJ(l) 


PARJ(2) 


PARJ(3) 


PAR J (4) 


PARJ(45) 


mean value 
interval 


0.091 
± 0.008 


0.31 
± 0.04 


0.35 
± 0.08 


0.045 
± 0.012 


0.50 
± 0.50 



Table 9: Parameter ranges used for random generation of Monte Carlo parameter sets without 
the popcorn effect. 

and an antibaryon share two valence quark-antiquark pairs. More generally, one expects 
any fragmentation model with very strong valence quark correlations between baryons and 
antibaryons to be disfavored. 

Appendix: Tuning of the Monte Carlo generator 

The tuning started with the PYTHIA parameters of [24], listed as "standard OPAL tune" in 
table 8. In a first step, the parameters, except the popcorn parameter, were varied individually 
and minimum values of \ 2 were searched for. This process was iterated. Search ranges for all 
parameters were defined either by requiring a maximum increase in x 2 of 20 or, in the case of 
smaller changes, by allowing a parameter shift of ±100%. Finally, all parameters were varied 
randomly within these ranges to search for x 2 values lower than that of the solution already 
found. Between 200 and 300 random parameter sets were generated at fixed PARJ(5) for the 
final search and 10 5 events were generated per parameter set. The results for the selected 
parameters sets MCI to MC5 are based on 10 6 events. The contributions of the four classes 
of observables to \ 2 a t the minimum value 56 are about 40 from the baryon rates, 10 from 
the meson rates and 6 from the remaining six observables. 

The measured and simulated baryon rates per Z° decay are given in table 7, while the 
corresponding PYTHIA parameters are listed in table 8. The baryon was included in 
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the optimization to avoid its almost complete suppression. The results for the parameter sets 
MC3 and MC5 are similar to those in our previous publication [12], where the case without the 
popcorn effect was not investigated. The reproduction of the meson sector does not change 
much during the tuning and is therefore omitted from table 7. 

Finally, table 9 gives the parameter space for the study of the parameter dependence of 
F-jj. All cases with \ 2 < HO were kept, allowing average discrepancies between the data and 
the model as large as those obtained with the modified popcorn model. 
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Abstract 

Data collected around y/s =91 GeV by the OPAL experiment at the LEP e + e _ collider are 
used to study the mechanism of baryon formation. As the signature, the fraction of S~ hy- 
perons whose baryon number is compensated by the production of a £~, A or H~ antihyperon 
is determined. The method relies entirely on quantum number correlations of the baryons, 
and not rapidity correlations, making it more model independent than previous studies. The 
diquark fragmentation model without the popcorn mechanism is strongly disfavored with a 
significance of 3.8 standard deviations including systematic uncertainties. It is shown that 
previous studies of the popcorn mechanism are not conclusive if parameter uncertainties are 
considered. 
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1 Introduction 



The formation of baryons within a jet of hadrons has proved difficult to model and is still not 
well understood. While the shape of momentum spectra can be derived from QCD with the 
modified leading logarithmic approximation together with parton-hadron duality [1], more 
complex observables like correlations have not been derived from first principles. 

Several physical models like the thermodynamic [2], cluster fragmentation [3] or string 
fragmentation [4, 5] models, have been developed to describe baryon production. Of these, 
the most successful is string fragmentation, based on the creation of diquark-antidiquark pairs 
from the vacuum as illustrated in figure 1. In a production chain of particles along the string, 
a baryon and an antibaryon can be produced in immediate succession (figs. la,b) or else one 
or more intermediate mesons can be produced between them (figs. lc-e). This production of 
intermediate mesons, referred to as the popcorn effect [6, 7], is included as an option in the 
Monte Carlo event generators JETSET and PYTHIA [8, 9] and can be steered with a free 
parameter. 

Past experimental investigations of the popcorn effect made use of rapidity ordering of 
hadrons in the fragmentation chain. Intermediate mesons modify the rapidity difference be- 
tween associated baryons and antibaryons. Rapidity correlations between A A pairs produced 
in Z° decays have been studied by several LEP collaborations and the conclusion was that best 
agreement between the experiments and the JETSET Monte Carlo model [8] was obtained 
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Figure 1: Baryon production in the diquark model without (a,b) and with (c-e) the popcorn 
effect. 

with the popcorn effect included [10]- [12]. 

A contradictory result was reported by the DELPHI collaboration [13]. Their measure- 
ment is based on triple correlations between a proton, an antiproton and a charged pion close 
in rapidity. Because the popcorn effect enhances the pion density in the rapidity interval 
between the proton and antiproton, the minimum rapidity difference between a pion and a 
proton was compared for the particle orderings p7rp and 7rpp. The measurement indicated 
the absence of the popcorn effect. Insufficient modeling of the fragmentation dynamics could 
not, however, be excluded [13]. 

In this analysis, the popcorn mechanism is investigated in a different way, by tagging 
rare baryons and measuring the quantum numbers of correlated antibaryons. By not relying 
on rapidity differences, we obtain results that are more model-independent than previous 
studies. The data were collected with the OPAL experiment at the LEP e + e~ collider at 
CERN. Especially suitable for our purposes is the S~ hyperon. If a £~ is produced, its 
baryon number and strangeness are compensated either by an antinucleon and a kaon, as 
illustrated in figure la, or by an associated antihyperon. The case of A or S° production 
without the popcorn effect is shown in figure lb. 
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The same graph, with the bottom-most u-quark replaced by a d- or s-quark, describes 
associated S _ and H~ production. Associated production of a S+ and H° antihyperon with a 
£~ can only occur through the popcorn mechanism (figs. ld,e). This makes and 
correlations ideal tools to study the popcorn effect. Unfortunately, the production rate of E~ 
hyperons and the probabilities for simultaneous reconstruction of S~ and £+ or £~ and H° 
particles are too small to make such an analysis feasable, given the available data statistics. 
The equivalent analysis cannot be performed with tagged A hyperons because in this case 
£+ and H° antihyperons can be produced without the popcorn effect. Alternatively one can 
measure the fraction 

of £~ hyperons accompanied by a E~ (F E _|p), a A (i 7 ^- j) or a H~ (F s _ — ). These 
correlations can occur in the popcorn model but are more likely in the diquark model (fig. lb) 
and their rate is thus a sensitive measure of the baryon production mechanism. 

The exact definition of the three correlations F s _ j in the sum (1) needs to account for the 
possibility that an event may contain more than one £~ hyperon, other additional hyperons, 
or more than one antihyperon. More generally, one can consider an arbitrary particle k and 
its antiparticle k. If the number of £~ hyperons in an event is larger than 1, all combinations 
of £~ k and E~ k pairs are counted. Denoting the total rates of £~-antiparticle and £ _ - 
particle pairs by i? s _ ^ and i?£-,fc> respectively, and the total S~ rate by R^-, the fraction 
F s _ £ can be written as 

This definition implies that the antiparticles k in an event enter n-^n^- times. The X~ Yr 
pairs are counted n s - (n s - — 1) times. In the data analysis, this multiple counting is not 
an issue because the number of reconstructed E~ hyperons per event is very small. Baryon 
number conservation ensures that Ylk^-k = ^> ^ ^he sum ex tends over all antibaryons, 
including the antihyperons and antinucleons. It is the understanding throughout this paper 
that the charge conjugated channels are included. Thus, the rates R^,-,k contain all like-sign 
pairs X~ k and £~ k, the rates ^ contain all unlike-sign pairs E~ k and S~ k, and £~ 
antihyperons are included in . 

It is evident from figure lc that kaons created through the popcorn effect reduce Fjj. 
The PYTHIA Monte Carlo program without the popcorn effect, tuned to reproduce the 
observed baryon rates and momentum spectra, predicts Fjj 0.9, as will be shown later. 
The other extreme is a model in which baryon number and strangeness are compensated 
statistically, i.e. from conservation laws alone. The ratio of weakly decaying hyperon to total 
baryon production in Z° decays is approximately 0.22 [14]. Because there must be at least one 
antibaryon in the rest of the event, if the T,~ is detected, and the tagging biases the number of 
strange valence antiquarks, the lower bound will be a bit larger: F-jj > 0.22. A more rigorous 
calculation can be performed with the thermodynamic model for particle production. The 
advanced version of the model, based on the microcanonical ensemble, yields F-jj = 0.23 for 
an initial system without strangeness [15]. 
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2 Event selection 



2.1 Event topologies 

In any experiment with a sufficiently large tracking device S~ hyperons can be identified 
by track kinks from X~ — > n-zr" decays. At LEP energies, the efficiency is much less than 
100%, because the decay vertex lies often outside the fiducial volume for its reconstruction. 
In principle, correlated decays £~ — ► n-7r + and H~ — > A7r + could be reconstructed using track 
kinks, too, but the efficiency for the exclusive reconstruction of £ _ S _ or £~ H~ pairs is 
small. 

In this work, only the S~ hyperons were reconstructed exclusively. Two signatures for 
correlated antihyperons were used: 

1. A hyperons were reconstructed by analyzing their so-called V° topology from the decay 
A — > p7r + . The decay vertex and the flight direction of the A candidates allow the 
A impact parameters do with respect to the beam line to be computed. Direct A 
production, including the contribution of decays from E° hyperons, is characterized by 
low impact parameters do, while large impact parameters indicate a preceding weak 
decay and are a signature for S decays. 

2. Charged pions with significant impact parameters are a signature for weak decays of 
arbitrary antihyperons. An inclusive sample of tracks with large do values, consistent 
with a pion interpretation, was selected. Throughout this paper, this data set is referred 
to as the sample of displaced tracks. This sample has a large background, and different 
antihyperon species contribute with different weights, because the number of decay 
pions per antihyperon is 1 for £~, 0.64 for A and 1.64 for chain decays of S - . 

The correlated E~A candidate sample gives the numbers of true S _ A and £~ H~ pairs. 
The fraction F^r can be extracted from the displaced track sample with a subtraction pro- 
cedure, using this information as a constraint. 

2.2 Experiment and data sets 

All data taken by the OPAL experiment [16] in the Z° energy region during the years 1991 to 
2000 were analyzed to measure the correlation. The OPAL experiment had nearly complete 
solid angle coverage and excellent hermeticity. The innermost part of the central tracking 
detector was a high-resolution silicon microvertex detector, which immediately surrounded 
the beam-pipe [17]. It was followed by a high-precision vertex drift chamber, a large- volume 
jet chamber [18], and z-chambers, all in a uniform 0.435 T axial magnetic field. In this work, 
the outer detector parts as well as the forward detector system were needed for triggering 
and identification of multihadronic events only. The criteria for multihadronic event selection 
have been described elsewhere [19]. 
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The present analysis is entirely based on the central tracking system [16, 18]. For candi- 
dates to be accepted, all central wire chambers and the microvertex detector were required 
to be fully operational. The data sample for this analysis consists of 3.97 million events. 

The identification of strange particles is based on earlier work [20, 21]. The innermost 
sense wires of the jet chamber had a distance of 25.5 cm from the beam spot, the wire-to- wire 
distance was 1 cm and there were 159 sensitive radial layers. The requirement of the pattern 
recognition program was the existence of at least 12 hits in the jet chamber. This makes it 
possible to identify S~ hyperons with decay lengths larger than 36 cm. 

The z coordinate, along an axis parallel to the electron beam, was measured with a 
precision of 700 fim with the stereo wires of the vertex chamber and 100 to 350 fim with the 
z chambers. For tracks leaving the drift chamber at the side cones, the z coordinates of the 
exit points can be computed from the radius of the last wire with a hit. At the beam spot 
a constraint can be set using the bunch length of the beam. Inside the jet chamber, the z 
coordinates were measured with the charge division method with a resolution of 6 cm. This 
is one limiting factor for the kinematical reconstruction of S~ hyperons to be discussed later. 

The quality of the impact parameter measurement is directly connected to the detector 
resolutions in the (r, ip) plane perpendicular to the beam axis. These resolutions were 5 fim 
to 10 fim for the microvertex detector [17], 55 fim for the vertex chamber and, on average, 
135 fim for the jet chamber [16]. 

To study the detector response, the selection was applied also to Monte Carlo samples 
used before at OPAL. These were generated with the JETSET7.3 and JETSET 7.4 programs, 
followed by a full detector simulation [22]. The steering parameters for the generator are 
given in [23, 24]. The subset of parameters relevant for this analysis is discussed in section 4.2 
and the numerical values are given in table 8 (appendix). In total, the Monte Carlo samples 
consisted of 4.65 million multihadronic Z° decays. 

2.3 £ Selection 

The production cross section of S~ hyperons has previously been measured by OPAL [20]. 
Here the selection criteria to find track pairs forming a decay vertex were slightly changed 
relative to [20] to improve the sensitivity of the analysis to correlations. The applied cuts 
are summarized in table 1. They accommodate pattern recognition tolerances and define a 
fiducial jet chamber volume to guarantee minimal hit numbers for both tracks and to remove 
background from the end plates. The ratio h/pi in the last line of table 1 is proportional to 
the decay time in the rest frame of the decaying particle. High values are rejected to reduce 
the substantial background from kaon decays. 
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variable 


condition 


transverse momenta of the decaying 
and the daughter particle 


Pti > 0.15 GeV/c, pT2 > 0.10 GeV/c 


total momentum of the decaying particle 


pi > 2 GeV/c 


impact parameter of the primary track 


do < 5 cm 


endpoint radius of the primary track 


te\ < 170 cm 


maximum gap between the primary and 
the secondary track in the (r, <p) plane 


d\2 < 10 cm 


number of secondary tracks 


1 


particle charges 


equal 


arc distance between the intersection point 
and the end points of tracks 1,2 in (r, ip) 


fii > —7 cm, fi2 < +7 cm 


radius of the intersection point 


fvtx > 35 cm 


distance of the intersection point 
from the end plate in z direction 


Azytx > 40 cm 


X 2 for agreement of the two tracks in the 
(r, ip) plane, using (3 x 3) error matrix 


X 2 > 2000 


ratio of the track length to the 
momentum of the decaying particle 


h/ Pl < 15 cm/(GeV/c) 



Table 1: Selection of track pairs forming a decay vertex (for a more detailed description of 
the variables, see [20]). 

Monte Carlo studies show that the following processes have to be considered as sources 
for S~ candidates, others being negligible [20]: 

1. YT -► nvr'; 

2. S+ ^n^-; 

3. S - -► Avr"; 

4. K~ — > 7r~+neutrals; K~ — > /i~+neutrals or K~ — > e~+neutrals; 

5. secondary reactions in the detector material and fake background. The last sample 
consists mainly of scattered particles. 

To determine the relative contributions of these processes to the observed data sample, an 
unfolding procedure very similar to that used in [20] was used. The invariant mass mj- and 
the decay angle 6* of the pion candidate in the rest frame of the hypothetical £~ particle were 
computed, assuming that the unseen neutral particle is a neutron. The angle 9* is defined 
with respect to the flight direction of the £~ candidate at the decay vertex. 

Figure 2 shows scatter plots of these variables for the data and three Monte Carlo sources 
of events with kinked tracks, namely events with H~ or K in the final state. The 
background 5, not shown in figure 2, is a smooth function of both observables and is largest 
in the backward direction cos#* = — 1. 
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Five two-dimensional regions, denoted by a, ...e, were introduced in the (m^- , cos 9*) plane 
to enrich dedicated kink sources (see figure 2). The bin number of a £~ candidate was used 
as an observable in the unfolding procedure to be described in section 3. 

In total, 16790 candidates were found in the two-dimensional plane in the mass range 
from threshold to 1.5 GeV. In the Monte Carlo sample 18754 kinks were identified. 




Figure 2: Definition of two-dimensional bins to disentangle S~ hyperons from background. 
The reconstructed mass is plotted versus the pion emission angle in the hypothetical center- 
of-mass system. Top left: data. Other plots: Monte Carlo simulations for three particle 
classes as indicated. The labeling of the two-dimensional bins is shown in the bottom left 
plot. The preferential bins for the kink sources are: a+d for £~ and £+, b for E~ , b,c and d 
for K~ and e for background. 
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2.4 Selection of pions from weak decays 



Two pre-cuts were applied to select tracks with a good reconstruction quality: the transverse 
momentum with respect to the beam had to be larger than 0.15 GeV/c and the track angle 
at the beam spot relative to the beam direction was restricted to the region | cos6>| < 0.80. 

To remove charged particles from charm or bottom decays, the impact parameter do^ 
with respect to the primary vertex was required to be larger than 0.2 cm. This cut is the 
essential condition to define the sample of displaced tracks. A lower value would be sufficient 
but does not improve the accuracy of the correlation analysis. 

Pions were enriched using the specific energy loss measurement of the central drift cham- 
ber [26]. A weight w^E/dxi^) was defined as the probability that the energy loss dE/dx of a 
pion deviates from the median value < dE/dx > (tt) by more than the measured difference 
from the median value. The applied condition was w^E/dxi^) > 0.02 and the number of hits 
contributing to this measurement had to be at least 20. 

Finally, the angle between the momenta of the pion candidate and the £~ candidate at 
the primary vertex was required to be less than 90 degrees. This hemisphere cut is motivated 
by the fact that it rarely happens that baryon number and strangeness are compensated by an 
antihyperon in the opposite event hemisphere. The cut reduces the combinatorial background 
by a factor 2. 

In total, 9965 correlated like-sign £ - 7r~ and 11951 unlike-sign £~ tt + pair candidates 
were selected with these cuts. The corresponding results for the Monte Carlo sample are 
10769 and 13818, respectively. If the total number of Monte Carlo track kinks is scaled to 
the observation, the number of like-sign pairs in the data is well reproduced by the Monte 
Carlo, the difference being (+3. 4 ± 1.4)%. The observed unlike-sign minus like-sign difference, 
however, deviates from the prediction by (—28 ±7)%. This deficit, already visible at raw data 
level, indicates that the Monte Carlo sample contains too many correlated antihyperons and 
is the basis for the final result of this paper. 

2.5 A Selection 

The selection cuts to find A decays in the central drift chamber have been described in [21]. 
The preselection cuts were relaxed from those of [21]. In this analysis all candidates with 
reconstructed masses between the p7r threshold and 1.20 GeV/c 2 were accepted. This larger 
mass window was needed to study the non-A background. 

Two criteria were added. If the reconstructed A flight path points back to the (r, ip) 
position of the hypothetical £ decay kink within 2 degrees, it was assumed that the kink 
originated from a H - — > Air~ decay and the A candidate was dropped. The cut removed 
approximately 3/4 of the like sign H~ A (S - A)-pairs and reduced the self-correlation of H~ 
particles with their own decay A's accordingly. 

Secondly, the hemisphere cut applied to the decay pions was also applied to the A's. The 
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angle between the flight directions of S and A at the primary vertex was required to be less 
than 90 degrees. 

For the correlation analysis, the mass window was reduced to a ±10 MeV wide interval 
around the true A mass. Totals of 276 £~A (£~ A) and 604 X~ A (£ _ A) pair candidates 
passed all selection cuts. In the Monte Carlo sample 284 like-sign and 721 unlike-sign pairs 
were found. The observed unlike-sign minus like-sign difference is smaller than the Monte 
Carlo prediction by (16 ± 11)%, if normalized to the number of observed kinks. 



3 Correlation analysis 
3.1 Unfolding of kink sources 

From the Monte Carlo sample one gets, for every kink source i and every two-dimensional 
(m s -,cos^*) bin j, the number of accepted events K^ MC \j). The populations of the bins 
j =a,...e in figure 2 are sensitive to the invariant mass and 9* resolutions, which in turn depend 
on the z resolution of the central drift chamber. The z coordinate, however, is not well modeled 
in the Monte Carlo simulation. Because the track end points in the jet chamber are well known 
by other measurements, this mismodeling at the decay vertex gives the dominant contribution 
to the systematic error of the X - rate [20]. To correct the Monte Carlo program for it, the z 
components of all S momenta were modified according to p z , n ew = p z ,rec + c- (p z ,rec —Pz,true), 
where p z ,true is the true momentum from the MC generator and p z ^ec the reconstructed 
momentum. The constant c is one common factor to be determined in the analysis. 

The total number of kinks K(j) expected in bin j is given by 



where Ndata and Nmc are the total number of multihadronic data and Monte Carlo events, 
respectively. Incorrect Monte Carlo rates are corrected for by the five scaling factors £j. If 
they are known, the true production rates per multihadronic event, can be computed for 
all sources i, for example 



Contrary to [20] the factors £j were treated as momentum independent, because the analy- 
sis of [20] had shown that the modeling of the spectral shape was satisfactory. The contri- 
bution from £+ hyperons to K(j) is less than 20% of that of the S~ particles [20]. The 
ratio of the genuine production rates, approximately 1 due to isospin symmetry, was fixed 
to the Monte Carlo prediction, so that £^q7 = . Four scaling factors £j were thus left for 
adjustment. They were computed with a x 2 -fit of the measured kink rates to equation (3) 
for a given value of c. 




(3) 




(4) 



j 
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It was then checked whether the Monte Carlo simulation reproduces the reconstructed 
mass and cos 9* distributions. Differing from our previous analysis [20], the proportionality 
factor c was chosen to get the lowest x 2 -sum for both distributions. The best overall agreement 
was found for c = 1.35 with a one-sigma interval ranging from 1.32 to 1.45. The final fit result 
is shown in figs. 3 and 4. The agreement is excellent. The error of c is not included in the 
statistical errors of the £j. It is treated separately as a contribution to the systematic error. 
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Figure 3: Measured £~ mass spectrum. Points: data. Histograms: results of the fit, ordered 
according to the sources for track kinks. The plotted errors are the statistical errors of the 
data. 
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The X - + X~ production rate was found to be R%- = 0.073 db 0.004. According to [20] a 
systematic error of ±0.009 has to be added. Within the total error, this result is consistent 
with our published value = 0.083 ± 0.011. The difference is mainly due to the modified 
treatment of the z-resolution. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos(0 ) 



Figure 4: Cosine of the center-of-mass pion emission angle. Points: data. Histograms: results 
of the fit, ordered according to the sources for track kinks. The plotted errors are the statistical 
errors of the data. 
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3.2 Correlated particle sources 



The correlated rates of kink-V° or kink-displaced-track pairs were measured as a function 
of the bin number j and the impact parameter do of the correlated particle relative to the 
beam line. The variable do contains information on the lifetime of the parent particle of the 
correlated particle in case of a preceding weak decay, and is needed to disentangle correlated 
A's from S's decaying into A's. 

In the Monte Carlo simulation, nine sources have to be considered for the displaced tracks 
and correlated V°'s: 



1. E- ^n^+; 

2. E+ — ► n-7r~ or p7r°; 

3. S~ — > Air + — > p7r + 7r + ; 

4. H° -► Avr° ^pvr+vr ; 

5. A — > p7r + , including A antihyperons from S° decays, but without the contributions from 
E decays; 

6. charged particles from K + decays; 

7. charged particles or V° configurations from K° decays; 

8. non-K° background of the V° topology; 

9. displaced tracks from secondary interactions, mainly scattering. 

All sources, except for the eighth one, contribute to the sample of displaced tracks. For V°-like 
events only the sources 3,4,5,7 and 8 are relevant. The sample of displaced tracks is enriched 
in pions but contains also a small fraction of protons and leptons. These contaminations are 
included in the Monte Carlo rates and are classified according to the above scheme. 

The combination of five kink sources with nine correlated particle sources leads to a total 
number of 45 classes of particle pairs. In addition, one has to distinguish between like-sign 
and unlike-sign pairs. Equation (3) can be generalized to the like-sign pair rate D like (j,d ) 
and the difference of the unlike-sign pair rate D unhke (j,do) and D hke (j,do), 

D^(j, do) = £ ril ■ ^ ■ D^ ke \j, do) ; (5) 

i,k 



D unlike^ do) _ D like {jj ^ = £ ^«) . ^ . {r> (M C, unlike) ^ ^ _ rfMCfike) ^ ^ (g) 

source i 



Nmc 

l,k 



The indices at the Monte Carlo rates £)WC,like) jj{MC,unlike) g p ec jfy ^ c kj n k 
and the correlated particle source k, respectively. The parameters rjf^ and rff^ are 90 scaling 
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• 

kink source 


correlated particle 


vj- 


£+ 






A 


K+ 


£- 


0.33 


0.014 


0.15 


0.01 


0.28 


0.04 


£+ 


0.01 


0.32 


0.01 


0.14 


0.28 


0.024 




0.34 


0.026 


0.15 


0.029 


0.36 


0.45 


K 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


0.03 


0.50 



Table 2: Fractions F- k as predicted by the PYTHIA 6.1 generator with the parameters 
from ref [24]- 



factors. In principle, they have to be extracted with a combined fit of the j and do dependent 
distributions of kink-track pairs and kink-V° pairs to equations (5) and (6). 

The Monte Carlo generator predicts the generic correlations F^ C \ defined with equa- 
tion (2), generalized to arbitrary sources i and k. Table 2 gives the results for hyperon- 
antiparticle correlations at generator level for the most recent parameter set used by the 
OPAL experiment [24]. 

If the scaling factors -qfl are known, the experimental results for the correlations can be 
computed 

Here, the scaling factors & from the kink fit (section 3.1) enter. 

The factors nfl are only needed to parameterize the statistical errors for the fit using 
equation (6). With very few exceptions, they are close to one and no physical result is 
extracted from them. 



3.3 Evaluation of the correlation matrix 

In spite of the large number of scaling factors rffl , a reliable fit of the data can be obtained. 

This can be seen from the integrals of the differences (D^ C ' unhke \j , do) — D^ C ' hke \j,do)) 
over the variables m£-,cos#* and do, listed in tables 3 and 4. The entries in the tables 
are normalized to the total differences, using rj^ = 1, and thus show the relative importance 
of the terms. Many correlations are very small. If not given in the tables explicitly, the 
contribution of a source is less than 2%. The entries in table 3 for displaced tracks account 
for 96% of the total rate difference. The contributions to the V° in table 4 add up to 104%, 
the excess being compensated by a small amount of S + A correlations with the opposite sign. 

The scaling factors for the four largest correlations in tables 3 and 4, r/°2 rj^l T , r/°2 
and n^l -r-, are determined with the fit. Three of them are needed to compute F-jj. 
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kink source 


displaced track source 


E- 




A 


K+ 


background 


xr 


(27 ±1)% 


(16 ±1)% 


(10 ± 1)% 




(17 ±4)% 




(12 ±1)% 


(6 ± 1)% 


(5 ± 1)% 




K 








(3±1)% 


kink background 











Table 3: Largest contributions to the correlations between kinks and displaced tracks. 



kink source 


V u source 




A 


V° background 


XT 


(20 ± 2)% 


(43 ± 3)% 


(14 ±3.5)% 




(6 ±1.5)% 


(21 ± 2)% 



Table 4: Largest contributions to the correlations between kinks and V° candidates. 



The data statistics do not allow to fit more than four parameters. The remaining 77's were 
thus fixed by symmetry considerations, isospin invariance or Monte Carlo studies. Systematic 
errors were assigned to them, if necessary. In the following, a few examples are described, 
preferentially the correlations in tables 3,4. 

For symmetry reasons, one has 

( a ) ( a ) to\ 

There is no model independent prediction for F~_ This parameter was set to the original 
Monte Carlo result. It follows then from equation (7) that 

4 a -^ = fe-- (9) 

It has been checked for various acceptable Monte Carlo generator tunings, described in the 
next section, that this procedure is valid within 20%. 

Non-negligible parts of the correlation are introduced by the correlated particle back- 
grounds 8 and 9; they are given as sums over the kink sources in tables 3 and 4. The origin 
of this effect is charge conservation in the events. Since the displaced-track-background con- 
sists mainly of scattered particles, it reflects the original particle charges. The kink selection 
introduces a charge bias for the rest of the event, which is visible in the ensemble of remain- 
ing charged particles on a statistical basis. This is true not only for scattered particles, but 
also for asymmetric fake V° candidates. Kinks and correlated particles are assigned to each 
other by chance. The condition ijftrack = f]fyo = £i was introduced, assuming a correct 
modeling of the background sources 8 and 9. This is justified because the equivalent relation 
for the like-sign background rates was confirmed with a fit to the data with the parame- 
terizations rjfjrack = &(track and r)f^ = CiCvo , the results being (track = 1-01 ± 0.04 and 
(yo = 1.03 ± 0.07, respectively. 

The remaining entry in table 3 is the charged kaon-kaon correlation. Its contribution is 
small because kaon decay inside the jet chamber is unlikely due to the long kaon life time. Since 
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the charged and neutral kaon production rates are almost equal and the hyperon rates are 
much smaller than the kaon rates, the fraction F K - K + is close to 0.5 so that rf^_ K+ = £k~- 

No model independent prediction exists for F K - ^, which was determined with the Monte 
Carlo event sample. A correlation potentially dangerous for the fit is the S~K + correlation. 
The corresponding asymmetry in table 2 is small, but actually it is the difference of much 
larger components. Strangeness conservation requires 

F £-,£- + F E-,S+ + 2 • F S-,H^ + 2 • F Z-*P + F Z-,A + 2 - F *~,K+ ~ 1 • (10) 

The parameters — and F^_ ^ appear with the weight 2 because the H-particles carry 

two units of strangeness. The contribution from hyperons is negligible. The factor 2 in 
front of F^-^k+ accounts for the K° contribution, which cannot be measured. Monte Carlo 
simulations with different model parameters showed that relation (10) is fulfilled within 3%. 
Together with the condition for baryon number conservation, 

F T,-^F + F E-,S+ + F T.-^BF + F E-,H0 + F S-,A + F ^~ ,P + F Y.~ ,n ~ 1, i 11 ) 

equation (10) gives the relation 

F Z-,K+ ~ \ ■ (*E",p + *E",n " F^_^r ~ ^) , (12) 

which shows the presence of large compensating terms. In the analysis, F^- x+ was therefore 
computed with the sum rule (10), taking the small correlations F s _ ^ and F s _ =o from the 
Monte Carlo generator. Similarly, F s - K + can be constrained by the equivalent equation for 
the E~, 

F E~,^ + F H-,S+ + 2 • + 2 • F S-,H0 + F S-,A + 2 " F S~ ,K+ ~ 2 . (13) 

The correlations of the S + hyperons are related to those of the T,~ hyperon by the isospin 
symmetry, for instance F s+ ^ = F s _ ^. 

The above relations allow either to replace rj^ by the four factors to be fitted, or to fix it 
and its contribution to the right hand side of equation (6). A simultaneous binned x 2 fit, using 
equation (6), was performed for the do distributions of the correlated-track and A candidates 
in the five (m s - , cos 6>*)-regions. In the A case, the reconstructed mass was restricted to the 
narrow interval given in section 2.5. In parallel, a fit of the like sign pairs to equation (5) was 
performed to determine (track = vftrack/^ an< ^ three normalization factors for the hyperon-A, 
K~-hyperon and (S~ or S~) — (£ _ or H - ) correlations. The statistical errors of the pair 
rates, needed to compute x 2 > were computed with the equations (5) and (6); they depend 
on the result of the fit. The fit was therefore done iteratively, setting the unknown rj- factors 
to unity for the first iteration. The statistical errors had do dependent fluctuations due to 
the limited Monte Carlo statistics. These were reduced with a one-dimensional smoothing 
algorithm described in [25]. The asymmetry fit resulted in a x 2 value of 144 for 146 degrees 
of freedom. The coefficients obtained are correlated, the largest correlation coefficient being 

£3-0.7 between — and r/ a ' ) — . 

'e-e- 'e-s- 
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3.4 Experimental results and systematic errors 

The adjusted do distributions are shown in figs. 5 and 6. The histograms give the contributions 
of the correlated particle sources; all kink sources and two-dimensional bins j =a,...e are 
combined. The errors of the data points are statistical and the corresponding errors of the 
Monte Carlo histograms are not shown. Fig. 7 gives the rate differences between the unlike- 
sign and like-sign kink-V° pairs as a function of the reconstructed A mass. Both distributions 
are very well described. 
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Figure 5: Impact parameters of the displaced tracks, correlated to £~ candidates. The 
difference between unlike-sign and like-sign combinations is shown. Points: data. Histograms: 
results of the fit, ordered according to the sources for displaced tracks. The plotted errors are 
the statistical errors of the data. 
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Figure 6: Impact parameters of the A candidates, correlated to S~ candidates. The difference 
between unlike-sign and like-sign combinations is shown. Points: data. Histograms: results 
of the ht, ordered according to the V° sources. The plotted errors are the statistical errors of 
the data. 
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Figure 7: Reconstructed masses of the A candidates, correlated to £~ candidates. The 
difference between unlike-sign and like-sign combinations is shown. Points: data. Histograms: 
results of the fit, ordered according to the V° sources. The plotted errors are the statistical 
errors of the data. 
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±0.013 


number of hits per track 
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+n no4 
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+0 006 


cos 9 distribution of tracks 
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A detection efficiency 
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±0.002 


charge asymmetry 
of detection efficiencies 


±0.002 


±0.001 


±0.001 


±0.002 


uncertainty of 
popcorn channels 


±0.008 


±0.005 


±0.004 


±0.007 


uncertainty of 
correlation 


±0.003 


±0.001 


±0.001 


±0.003 


uncertainty of 
K~A correlation 


±0.001 


±0.001 


±0.003 


±0.002 


final result 


0.17 ±0.11 


0.057 ±0.056 


0.25 ±0.08 


0.48 ±0.10 



Table 5: Experimental results and systematic errors. 



The fit results are listed in the first row of table 5. The errors are statistical and include 
the data and Monte Carlo contributions. Internally, the fit gets the A part of the correlation 
essentially from the correlated A sample, which sets also bounds on F^_ The 
correlation is computed from the displaced-track sample as a difference. Neither F^_ nor 
^=rr differ from zero in a statistical way, in contrast to the overall sum Fjj. The smaller 
error of the sum is due to the strong anti-correlation between F^_-^= and F^_ — . The final 
results, including the systematic errors, are given in the last row of table 5. 

The systematic errors, listed in table 5, will be discussed in the following. 

Mismodeling of correlated particle momenta. The momentum spectra of the correlated 
pion and A candidates are not well reproduced by the Monte Carlo program. The simu- 
lated momenta have to be scaled downwards by as much as 20%. The effect exists both 
for the like sign and the unlike-sign track and V° candidates. As a consequence, the 
measured fractions F i ^ are systematically too small. To find a correction, the detection 
efficiencies for the correlated particles were extracted from the Monte Carlo sample and 
the shifted momentum spectrum was folded with the efficiency function. This leads to 
an upwards correction of the measured correlations. A global correction of (6 ±3)% was 
applied. This correction, together with its error, is much larger than the uncertainties 
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due to miscalibrations of the particle momenta or the impact parameters. Therefore, 
no additional errors were assigned to the momentum and do selection cuts. 

z-resolution at kink vertex. In the ratio of a correlated rate to the single S~ rate the 
overall detection efficiency for E~ hyperons cancels. However, an uncertainty of the 
number of observed S~ hyperons arises from the unfolding procedure. The dominating 
error source is the mismodeling of the z-resolution. The correction factor c, introduced 
above, was conservatively varied between 1.25 and 1.5. For the individual correlations 
F-£~ t the obtained shifts were found to be non-parabolic functions of c. The maximal 
shifts of the results F s _ ^ and F-jj were taken as uncertainties. 

dE/dx calibration. Any antiproton impurity in the displaced track sample reduces the 
asymmetry from direct or indirect A decays. For a pure pion sample and a perfect 
calibration, the frequency distribution of the weight w^E/dxi^)-, defined in section 2.4, 
should not depend on w / dxi^) • A superimposed peak at w ^e / dxi 71 ) = is expected 
due to non-pions. The shape of this peak does not perfectly agree with the Monte Carlo 
prediction. A systematic error was assigned to the corresponding mismodeling of the 
antiproton rejection efficiency. 

Hemisphere cut. In principle, the results are corrected automatically for the hemisphere 
cut. No significant correlations were observed in the dropped hemisphere. A systematic 
error would appear, if the fragmentation model were incorrect. A 3% error was assigned 
to all correlations, based on the error of the number of tracks in the omitted hemisphere. 

Other systematic errors The number of hits per track is well modeled and introduces an 
error of 1% for the number of displaced tracks. Also the acceptance cut for the 8 angle 
plays a minor role only. From the difference between the angular distributions of data 
and Monte Carlo events an efficiency error of 3% was estimated. Systematic errors due 
to the cuts for A selection were already discussed in [21]. An uncertainty of ±3.3% 
was taken from that paper as fully correlated error for the A event sample. Differences 
between the detection efficiencies for particles and antiparticles would lead to a 
spurious asymmetry, if both the kinks i and the correlated particles k are affected. 
The effect is proportional to Ae s - • Ae^. The efficiency differences were extracted 
from the observed particle and antiparticle rates and the upper limits for the spurious 
asymmetries in table 5 were obtained. Furthermore, there are uncertainties due to the 
model dependence of the correlations which had to be subtracted. One half of the Monte 
Carlo baryon antibaryon pairs are accompanied by a popcorn meson. The 77^ -factors 
for the popcorn specific correlations and were varied between and 

2 and the fit repeated. The K~A correlation was varied by 50% and an uncertainty of 
20% was assigned to the correlation, fixing the small residual contribution from 

E~ self correlations. As already mentioned, the definition of F-jj introduces a small 
amount of double counting, both in the data and in the Monte Carlo simulation. This 
problem was studied with a toy Monte Carlo program and the only effects found were 
negligible corrections to the statistical errors of the fit. 

The correlated kink-track and kink-V° pairs should be concentrated in the S~ enriched 
bins a and d in the (m s - , cos 8*) plane shown in figure 2. As a cross check, a simpler analysis 
was performed, where the analysis was restricted to region a. All background asymmetries 
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and also the E~A correlation were subtracted as predicted by the Monte Carlo program. The 
final result including reevaluated systematic errors, F-jj = 0.472 ± 0.155, is fully consistent 
with the main result and demonstrates the absence of anomalies in the (m s -,cos6**) plane. 

4 Comparison with fragmentation models 

Before comparing the result for F-jj to the predictions of models, we adjusted the models, 
incorporated into the PYTHIA Monte Carlo event generator, to describe a set of observables in 
Z° decays. The optimization began with the PYTHIA steering parameters given in [24]. These 
were slightly modified to reproduce the newest experimental information on the baryonic 
sector, including data on AA correlations. The tuned Monte Carlo was then used to predict 
F-jj. We also used the models to study p7rp correlations more extensively than previously, as 
described below. 

4.1 Observables 

The input observables for the tuning fall into four categories: 

1. Eight baryon multiplicities in multihadronic Z° decay: p, A ++ , A, S + + X°, S _ , 
S* + + and VL~ . The production rates were taken from the compilation of the 
particle data group [14] . 

2. Proton and A momentum spectra. Baryon spectra have approximately Gaussian shapes 
if parameterized in terms of the variable £ = ln(l/x p ) [1]. We use the mean values and 
variances of the proton and A spectra as the observables. Deviations from the normal 
distribution are known, the true maxima being somewhat higher than the result of the 
fit [27, 31]. However, the Gaussian fit is an easy way to compare different experiments. 

The tabulated data [27]-[31] were fitted to Gaussian functions in the interval 1.2 < 
£ < 4.2. The situation for the protons is not satisfactory, the fitted maxima of the £ 
distributions varying from 2.79 [28] to 3.08 [31]. The values £ peak = 2.80 ± 0.07 and 
a = 1.11 ± 0.06 were obtained with a combined fit using all LEP and SLD data. The 
errors were not taken from the fit but conservatively estimated from the systematic 
differences between the £ spectra of the experiments. It should be noted that both £, pea k 
and a depend on the fit range due to deviations of the ^-distribution from the Gaussian 
shape. The agreement between experiments is better for A production, leading to values 
i peak = 2.62 ± 0.04 and a = 1.21 ± 0.04. 

3. AA correlations. The correlation was parameterized by two observables, the first one 
being the rate excess of A A pairs over A A + A A pairs per event, N^ r . The other 

observable is the mean rapidity difference Ay^ between the A and A after subtraction 
of like-sign pairs. The number of correlated A A pairs per event was taken from [12]. The 
mean AA rapidity difference Ay was computed from the data of [12], as the truncated 
mean for Ay < 3.0. This cut was introduced to suppress the contribution of A's from 
the opposite event hemisphere. 
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4. Nine meson production rates per multihadronic Z° decay: n + + n , 7T°, K + + K , K^, 
p + + p~, p°, if, K*+(892) + K*-(892) and K*°(892). These rates, taken from [14], were 
included in the tuning to protect the meson generation against parameter modifications 
steering the baryonic sector. 

4.2 Model Parameters 

The popcorn mechanism is incorporated in PYTHIA in two ways [9]. The first simple ver- 
sion was originally introduced in JETSET. Baryon production is controlled by the following 
parameters: 

1. the suppression of diquark-antidiquark production relative to quark- antiquark produc- 
tion, P{qq)/P(q) = PARJ(l); 

2. the suppression of ss production relative to uu production, P(s)/P(u) = PARJ(2); a 
tuning of this parameter was necessary, because the strange meson rates had to be 
readjusted; 

3. a double ratio involving diquarks containing s quarks, (P(us)/P(ud))/(P(s)/P(u)) = 
PARJ(3); 

4. the suppression factor for spin 1 diquarks, (l/3)P(«di)/P(udo) = PARJ(4); 

5. the popcorn parameter, which determines the relative occurrences of the baryon-meson- 
antibaryon and baryon-antibaryon configurations, PAR J (5); 

6. an extra suppression for having an s s pair in a baryon-meson-antibaryon configuration, 
PARJ(6); 

7. an extra suppression for having a strange meson in a baryon-meson-antibaryon config- 
uration, PARJ(7); 

8. a parameter which enters the exponent of the Lund symmetric fragmentation function 
for diquarks, PARJ(45); this parameter has an impact on the rapidity difference in 
baryon-antibaryon correlations. 

In the advanced popcorn scheme, a universal equation for tunneling from the vac- 
uum is applied to the generation of new partons and an arbitrary number of mesons can be 
created between a baryon and an antibaryon. The tunneling probability is proportional to 
exp(— [3 q ■ Mj_), where [3 q is a flavor dependent model parameter and M± the transverse mass 
of the created object. Only the first two and the last parameter of the above list are used in 
this scheme. There are three new parameters, two of them related to the tunneling formula: 

9. the tunneling coefficient for u-quarks, (3 U = PARJ(8); 

10. 5f3 = f3 s -f3 u = PARJ(9) 

11. an extra suppression factor for spin 3/2 baryons = PARJ(18). 
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Table 6: Comparison with fragmentation models. Fit quality x 2 , AA correlations and the E~ 
antihyperon correlations as defined in the text. The errors of the simulation are smaller than 
the last digit shown. 

4.3 Simulation Results 

For a given PYTHIA parameter set, the measured observables from section 4.1 were compared 
with the simulation results and a \ 2 was computed. It cannot be expected from a fragmen- 
tation model that all its predictions are correct to better than a few percent. To reduce the 
contributions of very accurately measured observables, the errors to compute \ 2 were thus 
taken to be at least 2.5%. Low \ 2 values were searched for with the method described in 
the appendix, for fixed values of the popcorn parameter. Many tunes of the generator have 
almost the same quality. Some PYTHIA parameter sets and the predicted baryon production 
rates are given in the appendix. Table 6 shows the x 2 ' s > the AA correlations and the £ _ - 
antihyperon correlations, as a function of the popcorn parameter, where the simple popcorn 
model is denoted MCI to MC4. 

The variable \ 2 ls an indicator for the quality of the baryon modeling. For the simple 
popcorn model, the lowest \ 2 value found was 56 for the 23 observables. This means that 
the Monte Carlo generator describes, on average, the observables roughly within 2 times the 
experimental errors or 5%, whichever is larger. 

The most important result of the simulation is that the overall quality of the description 
of the observables does not depend strongly on the popcorn parameter in the simple popcorn 
model. The AA correlation parameters are always reproduced within two standard deviations, 
whether the popcorn effect is switched on or off. The Monte Carlo parameter space examined 
here is larger than that in earlier studies. In view of the overall uncertainty it is not possible 
to reach a definite conclusion about the the popcorn effect by using AA correlations. 

As shown in the bottom line of table 6, the measured fraction F-jj is consistent both with 
the predictions of the original popcorn model with a large popcorn parameter (MC4) and the 
optimized advanced popcorn model (MC5), but is smaller than the Monte Carlo prediction 
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for zero popcorn effect. To investigate the stability of this last simulation result, the relevant 
parameters of the PYTHIA generator, especially the diquark fragmentation function, were 
randomly varied as described in the appendix. In these studies, Fjj values less than 0.86 
were not obtained. The experimental result deviates from this lower limit by 3.8 standard 
deviations. 



The advanced popcorn model has fewer parameters available for tuning and provides a 
significantly worse description of data as seen from the larger x 2 value. The larger x 2 f° r this 
model arises to a large extent from two well known facts. The absolute number of A A pairs 
is too large and the distribution of the A A rapidity differences is too narrow in comparison 
to the observation [12]. On the other hand, the average AA rapidity difference in the simple 
popcorn model with a very large popcorn parameter is too large, so that a combination of 
the two models might possibly describe the AA correlation well. 

To complete the comparison of correlations with results from Monte Carlo generators, the 
DELPHI p7rp [13] correlation was also investigated. This was done for rapidity ordered p7fp 
and 7rpp or pp7r particle configurations inside event hemispheres. The selection cuts and the 
definition of the minimum rapidity gap Ay m i n between a selected pion and the next proton 
were taken from [13]. The discriminating variable of [13] is the ratio of intensities 

( Urnin) ]V(p7fp) + iV(7rpp + pp7r) ^ ^ 

at the rapidity difference Ay m i n . A strong dependence on the popcorn effect had been seen 
by [13] at large values of Ay m i n . For three bins in the range 0.625 < Ay m i n < 1, the observed 
distribution agreed with a subsample of Monte Carlo events without the popcorn mechanism 
and disagreed with a disjunct subsample, containing popcorn mesons, by more than five 
standard deviations, averaged over the three bins. 

These results could be reproduced with the simulations described here: The R distribu- 
tions obtained with the parameter set MCI, but without detector corrections, agreed with 
DELPHI'S observation within 1.9 standard deviations, while there was disagreement between 
the data and Monte Carlo study MC2 by 5.2 standard deviations, averaged over the same 
three Ay m i n bins. Variations of the fragmentation model had not been studied in [13]. An 
increase of the fragmentation parameter PARJ(45) to unity, fixing the other parameters of 
the simulation MC2, reduces the difference between the data and the model prediction to 2.4 
standard deviations. Furthermore, it was found that the advanced popcorn model with the 
parameter set MC5 results in almost the same function R(Ay m i n ) as the model MCI without 
the popcorn effect. 

The modification of the parameter PAR(45) in the simulation MC2 increased the x 2 value 
in table 6 from 59 to 70 due to a shift of Ay^ to a value below the observation, without 
degrading the description of the other observables of section 4.1. In summary, these results 
indicate a high sensitivity of the rapidity correlations to the fragmentation dynamics. 
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Table 7: Baryon rates per multihadronic event. The errors of the simulation are smaller than 
the last digit shown. 



5 Conclusions 

Our investigations indicate that the fragmentation models have not yet reached a state where 
they can quantitatively describe the AA correlations, the ppzr correlations and the S~ antihy- 
peron correlation F-jj simultaneously. Neither the AA nor the ppzr correlations can provide a 
clear conclusion about the popcorn effect because they can both be described to an acceptable 
level by models either with and without the popcorn mechanism. The PYTHIA generator 
without the popcorn effect can reproduce both rapidity correlations simultaneously within 
two standard deviations. Both correlations are also in acceptable agreement with model 
predictions including the popcorn effect. However, we were unable to find a variant of the 
popcorn model that could simultaneously describe both types of correlations. 

In this work, the mechanism of baryon formation was studied by counting S~ hyperons 
and correlated antihyperons from hadronic Z°decays. The result F-jj = 0.48 ± 0.10, based on 
data taken by the OPAL experiment at LEP, favors a large popcorn parameter in the simple 
popcorn model and is also consistent with the advanced popcorn model. The fragmentation 
dynamics play no role here because the rapidity is not used in the analysis. Correlated particle 
momenta play an indirect role only, because they influence the detection efficiencies, but the 
final result contains a correction and a systematic error for mismodeling. 

Trivial correlations between hyperons and antihyperons based on baryon number and 
strangeness conservation, as predicted by the thermodynamic model, always exist. Due to 
the limited data statistics and the insensitivity to the dynamics, our result gives only a weak 
indication for non-trivial quark correlations between baryon-antibaryon pairs. Statistical 
models, only constrained by conservation laws, differ at most by 2.6 standard deviations. 

The result for Fjj deviates from the lower limit of simulations without the popcorn ef- 
fect by 3.8 standard deviations and thus demonstrates the existence of the popcorn effect. 
Any fragmentation model with very strong valence quark correlations between baryons and 
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Table 8: Parameter sets used for the simulations in tables 6 and 7. 
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Table 9: Parameter ranges used for random generation of Monte Carlo parameter sets without 
the popcorn effect. 

antibaryons is disfavored. 



Appendix: Tuning of the Monte Carlo generator 

The tuning started with the PYTHIA parameters of [24], listed as "standard OPAL tune" in 
table 8. In a first step, the parameters, except the popcorn parameter, were varied individually 
and minimum values of x 2 were searched for. This process was iterated. Search ranges for all 
parameters were defined either by requiring a maximum increase in \ 2 of 20 or, in the case of 
smaller changes, by allowing a parameter shift of ±100%. Finally, all parameters were varied 
randomly within these ranges to search for x 2 values lower than that of the solution already 
found. Between 200 and 300 random parameter sets were generated at fixed PARJ(5) for the 
final search and 10 5 events were generated per parameter set. The results for the selected 
parameters sets MCI to MC5 are based on 10 6 events. The contributions of the four classes 
of observables to \ 2 & t the minimum value 56 are about 40 from the baryon rates, 10 from 
the meson rates and 6 from the remaining six observables. 

The measured and simulated baryon rates per Z° decay are given in table 7, while the 
corresponding PYTHIA parameters are listed in table 8. The £l~ baryon was included in 
the optimization to avoid its almost complete suppression. The results for the parameter sets 
MC3 and MC5 are similar to those in our previous publication [12], where the case without the 
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popcorn effect was not investigated. The reproduction of the meson sector does not change 
much during the tuning and is therefore omitted from table 7. 

Finally, table 9 gives the parameter space for the study of the parameter dependence of 
F-jj. All cases with x 2 < HO were kept, allowing average discrepancies between the data and 
the model as large as those obtained with the modified popcorn model. 
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